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Abstract 

The properties of MogSg nanowires were investigated with ab-initio calculations based on the 
density-functional theory. The molecules build weakly coupled one-dimensional chains, like MogSe6 
and Mo6S9_ x I x , and the crystals are strongly uniaxial in their mechanical and electronic properties. 
The calculated moduli of elasticity and resilience along the chain axis are c\\ = 320 GPa and 
Er = 0.53 GPa, respectively. The electronic band structure and optical conductivity indicate that 
the MoeSe crystals are good quasi-one-dimensional conductors. The frequency-dependent complex 
dielectric tensor e, calculated in the random-phase approximation, shows a strong Drude peak 
in en, i.e., for the electric field polarised parallel to the chain axis, and several peaks related to 
interband transitions. The electron energy loss spectrum is weakly anisotropic and has a strong 
peak at the plasma frequency hv p ~ 20 eV. The stability analysis shows that MogSg is metastable 
against the formation of the layered M0S2 ■ 

PACS numbers: 61.46. +w, 62.25.+g, 73.22.-f, 73.63.-b, 78.67.-n, 71.15.Mb 
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I. INTRODUCTION 



Quasi one-dimensional materials like nanowires or nanotubes are attracting considerable 
interest, in particular because of their potential applicability in nanoelectronic and nanome- 
chanical devices. Common to all of them is very strong uniaxiality in the electronic as well as 
mechanical properties, which is related to weak van der Waals coupling between individual 
molecular chains. They are very elastic, the strain at the (ultimate) tensile strength exceeds 
that of steel often by more than two orders of magnitude. As a consequence, they show 
extraordinary tensile strength along the long molecule axis - although the Young moduli 
are of the same order of magnitude. Perhaps even more challenging are their electronic 
properties. By tuning their bandgap or conductivity, one can build a series of very efficient 
and compact nanoscale electronic devices. In the past decade the research has largely con- 
centrated on carbon-based nanotubes which often showed promising properties. Inorganic 
molybdenum-chalcogenide nanowires represent, as we shall see, a possible complementary 
or even alternative material with similar mechanical and electronic properties, but they are 
easier to prepare in clean, single-stoichiometry form. 

As early as in 1980 Potel et al. reported on the synthesis of M 2 Mo 6 S 6 (M = K, Rb, Cs). 1 
These compounds all crystallize in the P63/7??. space group with the MogSg chains oriented 
parallel to the hexagonal c axis, located along the rhombohedral c edges and with the alkali 
metals in the (2/3, 1/3, 1/4) and (1/3,2/3,3/4) positions in the unit cell. 

Molybdenum-chalcogenide crystals intercalated with Li (Lia;Mo 6 Se 6 , < x < 8) or with 
other metals were also investigated in the past, e.g., as solid-state electrodes for secondary 
(rechargeable) lithium batteries, as one-dimensional conductors 2,3 and even as superconduct- 
ing materials. 4 In the second example Li or Pb were needed to enhance the conductivity of 
nanowires. It seems that the extremely interesting mechanical properties of molybdenum- 
chalcogenide crystals were overlooked at that time. More recently, F.J. Ribeiro et al. 5 
investigated the elastic properties and the band structure of Li 2 MoeSe6 by performing the 
ab-initio total energy calculations. They pointed out the extraordinary high theoretical 
value of the Young modulus, comparable to the modulus of carbon nanotubes. 

Other Mo-based nanowires, composed of Mo, chalcogens (S) and halogens (I) to form 
MoeScj-xIx (Mo-S-I) have been investigated recently. 6 They have a skeleton of M06 octahedra, 
each dressed with six anions (S or I), and bound together into chains by three anions (again 
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S or I). Specific to these materials is growing in bundles of identical chains, large Young 
moduli, very small shear moduli, and easily controlled electronic properties. 

In this paper we concentrate on the Mo 6 S 6 (Mo-S) nanowires without and with interca- 
lated K atoms. Using the ab-initio total-energy and band-structure calculations we predict 
and discuss the stability, elastic and electronic properties of crystals made of Mo 6 S 6 . 

II. STRUCTURE AND COMPUTATIONAL DETAILS 

The crystal structure of Mo-S is shown in figure 1. The hexagonal crystals comprise one- 
dimensional chains oriented along the hexagonal c axis. Each chain is formed of a stack of 
M03 triangles with staggered rotations and dressed with three anions (S). In order to clarify 
the stability, possible structures and properties of Mo-S nanowires we performed extensive 
density-functional theory (DFT) analysis using the WIEN2k code. 7 The simulations were 
performed on a mixed basis set of augmented plane waves plus local orbitals (APW+lo) 
for low orbital momenta (£ < 2) and linearised augmented plane waves (LAPW) for all the 
higher orbital momenta. 8 The exchange and correlation potential was treated in the Perdew, 
Burke and Ernzerhof generalised-gradient approximation (GGA). 9 The muffin-tin radii were 
set to R$jg = 1.06 A for Mo, Rf T = 1.27 A for S and R^ T = 1.32 A for K, the kinetic 
energy cutoff was E^ ax = 12.3 Ry and the plane-wave expansion cutoff was E^ x = 196 Ry. 
The total energy was calculated on a mesh of 150 (9 x 9 x 18) Appoints in the irreducible part 
of the Brillouin zone. The problem of the above DFT method is inappropriate treatment 
of the weak long-range van der Waals interaction which is responsible for bonds between 
individual chains in such molecular crystals. As a consequence, the DFT treatment of 
interaction between individual chains is not very reliable. On the other hand, the analysis 
of K 2 Mo 6 S6 chains reveals that the calculated hexagonal lateral lattice constant a (at T = 
K, as usual in the DFT calculations) differs from the experimental one (at room T) only 
by 2.6%, see Table I. In any case, it is clear from both theory and experiment, that the 
lateral interaction between the chains is very weak. This is another reason why we shall 
concentrate on the intrachain structure and properties of single chains. 

In the next stage, the total energy was minimised with respect to the internal atomic 
coordinates and to the hexagonal lattice constants a and c by assuming that MogSg belongs 
to the same space group as K 2 MoeS6, i.e., P63/771 (176). 1 This assumption is corroborated 



3 



by weak coupling between individual molecular chains and on extensive analysis of Mo-S-I 
nanowires. 6 The space group PQ^/m allows rotational relaxation of molecular chains around 
the c axis but does not allow mutual translations or deformations of chains along the c 
axis. The results of the energy minimizations of MogSg and K^MogSg are compared with the 
available experimental lattice constants in Table I. The structure of MogSg is shown in figure 
1 and the interatomic distances in Table II. In K 2 Mo6S6 the calculated lattice constant a ca i c 
is 2.6 % larger than a exp , most probably because the DFT GGA overestimates the dispersive 
forces, acting between the chains. The nearest neighbouring distance between K atoms in 
K 2 Mo 6 S 6 (cInn = c) is close to the interatomic distances in metallic K (d — 4.53 A). 

The interatomic distances between the atoms of the same MogSg chain are up to 8% longer 
than the intrachain distances reported for M06S3I6 chains. 6 The in-plane Mo-Mo distance 
agrees, within the precision of our DFT calculation, with the experimental value for the 
room-temperature bulk metallic Mo (2.73 A). In the Supplement of Ref. 6 an agreement 
within 0.1 A between the room-temperature experimental distances in Mo 6 S3l6 and the 
T = distances calculated by the DFT was reported. 

III. RESULTS AND DISCUSSION 

To investigate the stability of the M06S6 and K2M06S6 nanowires we performed also the 
DFT simulation of the layered M0S2 compound and of atomic S and K in vacuum (i.e., one 
S or K atom in the unit cell 10.1 x 10.6 x 11.1 A 3 with periodic boundary conditions). The 
comparison shows that for Mo 6 S 6 the energy difference with respect to the more stable MoS 2 
is AE = £(Mo 6 S 6 ) -6£(MoS 2 ) +6E(S) = 28 eV (per one Mo 6 S 6 unit) and tells us that the 
Mo-S nanowires with Mo in oxidation state 2 are, like other Mo-S-I nanowires, metastable 
against MoS 2 with the Mo oxidation state 4. The stability of MogSe is comparable to the 
stability of already synthesised material M06S3I6 for which the calculated energy difference 
is AE = £(Mo 6 S 3 I 6 ) - 6£(MoS 2 ) + 9E{S) - 6E{1) = 32 eV. Per one atom, however, the 
energy difference is 2.33 eV for Mo 6 S 6 and 2.13 eV for Mo 6 S 3 I 6 . Therefore we expect Mo 6 S 6 
to be stable in the air, in analogy with the Mo-S-I nanowires. The already synthesized 
K 2 MoeS6 is stable in the air, 1 the corresponding energy difference to MoS 2 is 21.0 eV and 
the Mo oxidation state 2.33. 
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A. Elastic properties 



The Mo-based nanowires and nanotubes have large elastic moduli in the direction of the 
wires and very small shear moduli in general. 10 To investigate the mechanical properties of 
the Mo-S nanowires, we performed DFT simulations for a series of lattice constants c. Due 
to weak interchain coupling the effect of a on the total energy is very small and was neglected 
unless otherwise stated. The lattice constant a was thus kept fixed during expansion, but 
the atomic positions were fully relaxed. For larger lattice constants, when the nanowires 
start to break, we simulated double cells, composed of two unit cells in the c direction, so 
that each simulated cell contained 12 Mo and 12 S atoms and the space group was reduced 
to PI. In this way we reduced the effect of the interaction between breaking segments in 
a simulation with periodic boundary conditions. The results, presented in figure 2 show 
the behaviour of the stress an (cr n oc dE/dc at constant a) under tensile strain en- The 
initial slope yields for the elastic coefficient in the direction of the nanowire Cn ~ 320 GPa, 
which is larger than the Young modulus of steel but about three times smaller than for 
(single-wall) carbon nanotubes, see Table III. Upon increasing strain, the slope decreases 
and eventually changes sign. Usually we say that the system breaks after an reaches the 
maximum stress. Mo-S, however, behaves in a very different way. an decreases after the 
first maximum, but then it starts to rise again until it reaches a local maximum. Breaking 
of Mo-S nanowires is thus a two-stage process, see fig. 3. In the first stage, the Mo-Mo 
bonds at some place are stretched until they break. In the second stage, the nearby S atoms 
move into the empty space between the chains and bind together the two ends of the wires 
until, eventually, also these Mo-S-Mo bonds and thus the whole chain breaks at the second 
maximum in figure 2. Notice that the nonlinear regime below e n ~ 0.6 in fig. 2 is the 
consequence of the anharmonic contributions to the interatomic potential and is not related 
to plastic deformation which is absent in our DFT simulations. 

The modulus of resilience, i.e., the maximal energy that can be reversibly stored in an 
elastic medium, is in our case equal to 



The Mo-S nanowires have a very high modulus of resilience, E R ps 0.53 GPa (see Table 
III). This value is more than two orders of magnitude higher than for steel and already 
indicates possible applications of these nanowires: nanosprings. The ideal tensile strength, 
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i.e., the maximal stress is T rs 11 GPa, it is extremely high for an inorganic nanowire, but 
is by a factor of 3 smaller than for carbon nanotubes. Stretching the wire along the c axis 
by 6.4% causes a 3.7% elongation of the nearest interplanar Mo-Mo bond length and only 
a 2.6% contraction of the in-plane Mo-Mo bond length, see Table II. The wire breaks at 
en = 6.6 %. 

These are of course theoretical values. The experimental values of the tensile strength 
will be smaller due to lattice imperfections. Due to the almost one-dimensional structure 
and consequently small shear modulus, broken chains will diminish the tensile strength of 
an entire bundle of chains. 

B. Electronic properties 

The basis of electronic and optical properties of MoeSe is the electron band structure, 
shown in figure 4. For clarity, let us assign the sub-bands below -1.4 eV to the valence band 
and all the higher electron states to the conduction band. The Fermi energy then lies in 
the conduction band. Some of the sub-bands are doubly degenerated so that the conduction 
band carries 20 electrons per unit cell in total and all the valence bands together another 
52 electrons. These numbers agree with the total number of valence electrons in the unit 
cell: there are 12 atoms in the unit cell, each carrying 6 valence electrons. Due to large 
number of equivalent atoms in the unit cell it is not possible to uniquely assign particular 
sub-bands to individual atoms. Some general observations can nevertheless be made. The 
band at -12.5 eV (not shown in figure 4) is clearly the 3s— band of S atoms. Between -5 
and eV we have a mixture of Mo-Ad and S-3p bands, between and 5 eV the bands have 
a clear Mo-4d character, whereas above 5 eV the S-3d states start to prevail. In the whole 
energy interval there is also a substantial contribution of interstitial electrons to the total 
density of states. For comparison, the band structure of K 2 MoeS6 is shown in fig. 5. The 
intercalation of K atoms does not alter the general shape and anisotropy of the bands which 
is still governed by the intrachain wavefunctions overlap. K atoms add two electrons per 
unit cell which occupy one band and shift the Fermi energy up by about 1.5 eV. The above 
band structures are also similar to the calculated band structures of Mo 6 Se6 and Li 2 Mo6Se6 
nanowires. 5 

Dispersion of individual occupied sub-bands in the plane perpendicular to the wires is 
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very small, of the order 0.2 eV and is similar in magnitude to the dispersion in Mo-Se and 
other Mo-S-I nanowires. 5 ' 6 The reason for small dispersion lies in weak, mainly Van der 
Waals coupling between the wires. Interestingly enough, even intercalation with K does 
not alter the dispersion in this plane noticeably. The conduction electrons are localised 
on individual wires and the hopping rates between the neighbouring wires are small. The 
dispersion in the direction parallel to the wires is more system-specific. In MoeSe there are 
four sub-bands crossing the Fermi energy Ep, thus building four Fermi surfaces, shown in 
figure 4(a). The first surface is a small, weakly anisotropic sphere close to the T point which 
makes very small contribution to the charge density at Ep and to the conductivity. The 
main contribution to the static charge carrier transport and to static dielectric properties 
comes from electrons at the upper two Fermi surfaces which cross the r — A line close to 
the A point and are very flat as a consequence of anisotropic coupling between the wires. 
Thus, electrons in two sub-bands contribute to the static charge transport, the other either 
occupy sub-bands that don't cross the Fermi energy or have high effective mass. These two 
sub-bands have the dispersion in the direction parallel to the wires typical for nearly-free 
electrons. We find one electron sub-band starting at E(T) = —0.73 eV with an effective mass 
mjj 1 ~ 1.2m (m is the free-electron mass) and one hole sub-band starting at E(F) = 0.37 
eV with an effective mass mj| ~ 0.8m. The slopes at Ep of these two sub-bands give the 
electron group velocity ~ 0.4 x 10 6 m/s and the hole group velocity ~ 0.6 x 10 6 
m/s. The spin-orbit coupling - if added - neither splits the two subbands that cross the 
Fermi energy close to A nor does it alter noticeably the effective masses or the velocities 
of the two sub-bands. The values of m\\ and v\\ are of the orders of magnitude typical 
for metals, the Mo-S nanowires are thus candidates for one-dimensional ballistic quantum 
wires where the resistance is determined by the contacts of the wire. The ballistic regime is 
limited by the ability to produce defect-free nanowires, i.e., the length of the wire must not 
exceed the scattering length of the electrons. For a lifetime broadening of 0.01 eV (lifetime 
r ~ 7 x 10~ 14 s) this condition limits the length of the wires to ~ 40 nm. As one-dimensional 
conductors, the Mo-S nanowires are also candidates for a Peierls transition to an insulator at 
low temperatures. However, we cannot give an estimate here because we did not investigate 
the electron-phonon interactions. 

Let us compare the electron bands in Mo 6 S6 with those in MoeSg (three additional S 
atoms are ordered in triangles in the bridging planes between successive Mo 6 octahedra, see 
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Ref. 6) and in Mo 6 S 6 I 3 (in this case the Mo octahedra are "intercalated" with 3 I atoms 
instead of 3 S). In the latter two cases the sub-band dispersion in the z direction is smaller 
and the band structures have narrow gaps close but above the Fermi energy. The Mo-id 
electrons in MoqSq build conduction channels along the chains which are interrupted in 
case of Mo 6 S 9 and Mo 6 S 6 I 3 . As a consequence, the Mo 6 S 6 is expected to be a much better 
conductor than Mo-S-I chains of the type MoeX 9 (X is a combination of S and I). 

C. Optical properties 

During synthesis of Mo-S-I nanowires several different stoichiometries are obtained which 
then have to be purified in order to obtain a single-stoichiometry compound. Therefore 
it is essential to distinguish between them. To facilitate the characterisation of different 
stoichiometries we predict here also the optical properties of MoeSe- 

We start with the frequency dependent complex dielectric function, figure 6. The imag- 
inary part of the dielectric tensor e in the limit q — > is calculated with the WIEN2k 
code on a tetrahedral mesh of 300 A;— points in the irreducible part of the Brillouin zone 
and in the random-phase approximation. 13 The real part is then obtained from Im(£) with 
the Kramers-Kronig relations. In our case the e tensor is diagonal with the components 
£zz = £\\ and e xx = e yy = e±. The dielectric tensor has a Drude peak, associated with the 
transitions between the conduction electrons close to Ep, which is more pronounced in e\\. 
The anisotropy in e(E — > 0) is a direct consequence of the shape of Fermi surfaces. Above 
the Drude peak, has two pronounced absorption peaks at 3 and 5 eV whereas e\\ has a 
structured absorption peak between 2 and 2.8 eV. These peaks are related predominantly 
to interband transitions between different S— 3p (below Ep) and Mo— 4c? (above Ep) sub- 
bands. Above 5 eV the transitions are mainly to/from the S— 3d sub-bands. In all the cases, 
however, the influence of interstitial electrons is large. Coming from the high-frequency side, 
Re(e) changes sign (vanishes) at the plasma frequencies hio v ^ = 20.1 eV and ftw p \\ = 20.0 eV. 
The plasmons are longitudinal collective excitations of all valence and conduction electrons. 
u p agrees well with the estimate oj^ = ne 2 /e m (e is the electron charge, m its mass and e 
the permittivity of free space). With n = 72 valence electrons per unit cell of the volume 
V = 317.6 A 3 we get ftw p = 17.6 eV. The plasma frequency can be also tested with the 
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/-sum rule, 

71 f°° 

= J Q hn{eii(w)]wdu- (2) 

Numerical integration in the range from to 41 eV gives hw± = 18.5 eV and fkv\\ = 18.9 eV. 
At the plasma frequency the effect of the anisotropy in the band structure on the dielectric 
tensor is small and e is almost isotropic. 

From e we get the optical conductivity Re[cr(a>)] = e ulm[e(u)], see figure 7. The static 
{uj — > 0) peak is lifetime broadened due to the charge- carrier scattering mainly on impurities 
and imperfections in the crystal lattice. In the figure, the Drude peaks are broadened with 
T|| — r_L = 0.1 eV, which corresponds to a lifetime r pa 0.7 x 10 -14 s. With such F the 
static conductivity along the wires is cry ~ 2 x 10 4 S/cm and the predicted anisotropy in a 
is 100 : 1. At higher frequencies (fuv > 5 eV) the anisotropy in e gradually disappears since 
the electrons are excited also above the ionisation threshold of 4.36 eV into the unbound 
(free) electron states. 

The complex index of refraction (n + ik) with the components 

n u = y H+|Sj (3) 

and 

ku = yHZ|S (4) 

is shown in figure 8. It has similar trends as the dielectric function from which it is derived. 

Finally, the electron energy loss function Lu(u>) = — Im (l/eu(u>)) is shown in figure 9. 
The pronounced peak at ~ 20 eV lies very close to u p and is intimately related to the plasma 
excitations; 14 the effect of intraband transitions on this peak is small. 

IV. CONCLUSIONS 

The aim of the paper is to stimulate synthesis and research on Mo 6 S 6 and M 2 Mo 6 S 6 
(M = alkali metal) which - according to our theoretical predictions - are good quasi-one- 
dimensional conductors and at the same time have high moduli of elasticity and resilience 
along the wires. This combination of electronic and mechanical properties makes them 
unique in the area of one-dimensional nano materials. The chains made of triangular Mo 3 , 
surrounded with three S, build conduction channels along the hexagonal c axis with very 
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weak "leakage" of the current between individual chains. The intercalated alkali metals 
contribute additional charge carriers and thus cause a shift in the Fermi energy but do not 
build their own conduction channels. 

Synthesis of MogSg should follow the guidelines to make the MogSeg 2 whereas M 2 MoeS6 
has been synthesized quite some time ago together with M 2 Mo 6 Se 6 . 1 Synthesis might be 
difficult since in particular the compounds with S are metastable against the much more 
stable layered MoS 2 . Intercalation with alkali metals, if necessary to stabilise the compound, 
does not change the described properties qualitatively. 

Mo 6 S 6 and M 2 Mo 6 S 6 could be challenging alternatives to other nanomaterials like carbon 
nanotubes. The compounds have similar mechanical properties and thermal stability, they 
are stable to about 700 K in the air (the thermal stability of MogSg is believed to be 
comparable to the stability of other already existing Mo-S-I compounds). The advantage of 
Mo 6 S 6 and M 2 Mo 6 S 6 could be that they are easier to synthesise in a clean, one-stoichiometry 
form and that they have better conductivity than carbon nanotubes, due to their purely 
molybdenum metallic chains. Mo 6 S 6 is, according to the DFT results, also better conductor 
and has better mechanical properties than the already synthesised Mo-S-I nanowires like 
Mo 6 Sg_ x I x . One dimensional nature of the systems opens also the possibility of ballistic 
conduction or Peierls instability. 
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TABLE I: The calculated lattice constants (a, c) ca i c of MogSg and K2M06S6 and the experimental 
lattice constants (a, c) cxp of K2M06S6. 1 



Mo 6 S 6 K 2 Mo 6 S 6 


a calc 


9.2 ±0.1 A 9.05 ±0.1 A 


^exp 


8.82 A 


Ccalc 


4.35 ± 0.05 A 4.43 ± 0.05 A 


Cexp 


4.44 A 



TABLE II: Calculated ground-state interatomic distances in MoeSe- The in-plane distances are 
between the atoms of the same Mo triangle with the adjacent three S atoms and the distances 
between the planes are between atoms in nearest neighbouring planes. The numbers in the last 
column are for c = 4.63 A, (en = 6.4 %), i.e., close to the point where the chains start to break. 



Bond 




Distance 




In equilibrium Strained 


Mo-S In-plane 


2.49 A 


2.46 A 


Mo-S Between planes 


2.62 A 


2.65 A 


Mo-Mo In-plane 


2.74 A 


2.67 A 


Mo-Mo Between planes 


2.69 A 


2.79 A(2.782, 2.798) 


S-S Between planes 


3.60 A 


3.66 A 


S-S Closest interchain 


4.22 A 


4.29 A 
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TABLE III: Elastic properties of nanowires, carbon nanotubes (single-wall, SWCNT and multi- 
wall, MWCNT) and steel for comparison. Y is the Young modulus, Er modulus of resilience, and 
T the tensile strength. The data for uniaxial systems are along the wires or tubes. The ideal tensile 
strength of MogSee is overestimated because it was assumed that the crystal stretches uniformly. 



Y (GPa) E R (GPa) T (GPa) 
Mo 6 S 6 a 320 0.53 11 

K 2 Mo 6 S 6 a ~230 

Mo 6 Se 6 6 250 3.1 26 

SWCNT C 1000 5 to 10 ~30 
MWCNT < 1800 

Steel 210 ~3x 10~ 3 0.4 

"Calculated in this paper 

b Calculated using the data for an isolated molecular chain 5 
Compilation of data from Refs. 11 and 12 
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FIG. 1: The crystal of MogS6. (a) top view (along the [0001] direction) and (b) side view (along 
the [1210] direction). Grey (dark) spheres: Mo, yellow (light) spheres: S. The MogSg molecules 
form chains along the [0001] direction and a bundle of chains forms a nanowire. The separation 
between the chains is large and the interaction between them is weak, the MogSg crystal is strongly 
anisotropic. In K2M06S6 the M06S6 chains are rotated by ~ 20° in the anticlockwise direction and 
the K atoms are in the (2/3,1/3,1/4) and (1/3,2/3,3/4) positions. 
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FIG. 2: Elastic properties of MogSg nanowires. The elastic strain en is for a unit cell, doubled 
along the c axis. At the first and second maxima the Mo-Mo bonds and the Mo-S-Mo bonds break, 
respectively. The initial slope of the curve defines the elastic modulus en = 320 GPA, and the 
grey area the modulus of resilience Er = 0.53 GJ/m 3 . 




a b c 



FIG. 3: Stretching of MogSe nanowires is a two-stage process, as shown on this graph with two unit 
cells in the direction of the chain, (a) Fully elastic regime before the Mo-Mo bonds at a particular 
chain segment start to tear apart (c = 4.60 A), (b) The Mo- Mo bonds between two neighbouring 
unit cells start to break and the S atoms move into the empty space between them by keeping the 
two ends of the broken chain together via four Mo-S-Mo bonds (c = 4.975 A), (c) Eventually also 
the Mo-S-Mo bonds break (c = 5.29 A). 



15 



r E M K A T A A "0 20 40 

DOS 

b c 

FIG. 4: (a) The Fermi surfaces in the Brillouin zone of MoeSg are very anisotropic and reflect 
the fact that the coupling between the molecular chains is weak. As a consequence, the electron 
group velocity points mainly parallel to the c axis. The letters show the positions of symmetry 
points in the Brillouin zone, (b) Electron band structure along the main symmetry directions. The 
dispersion of occupied states in the lateral direction is small, reflecting weak coupling between the 
chains. The radii of the circles indicate the Mo— Ad character of the bands. The main contribution 
to the static charge-carrier transport comes from two bands, crossing the Fermi energy close to the 
A point, (c) The total electron density of states per eV in a unit cell. 
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FIG. 5: Electron band structure of K2M06S6 along the main symmetry directions. The main 
difference to the bandstructure of M06S6 is the shift in Ep. 
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FIG. 6: Energy dependence of the complex dielectric functions en and e±, broadened by a 
Lorentzian with T = 0.1 eV. The dielectric function en has a pronounced Drude peak, typical 
for intraband transitions in metals, and an absorption peak at 2.8 eV. e± shows two pronounced 
absorption peaks which are associated to the interband transitions and are closely related to the 
maxima in the joint density of states, which means that most of the sub-bands in the valence and 
conduction bands are a mixture of several symmetries, they have no unique band character. As a 
consequence, the selection rules are not very selective. 




E(eV) 



FIG. 7: Spectral dependence of the optical conductivity components Re(<7||) and Re(<7j_). Clearly 
seen is strong anisotropy in the static conductivity. The otherwise diverging peak at E = was 
damped with Tu = T± = 0.1 eV to mimic the charge-carrier scattering. 
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FIG. 8: Complex index of refraction (n + ik), consisting of ny and n± (upper panel) and extinction 
coefficients k\\ and k± (lower panel). 
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